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The compression of a low-density copper foam was simulated with a radiation, hydrodynamics 
code. In one simulation, the foam had a density of po = 1-3407 g/cm 3 , 15% the density of copper at 
standard temperature and pressure, and was composed of a tangle of standard density cylindrical 
copper filaments with a diameter of 4.0 x 10 -7 cm. In another simulation, the foam was a uniform 
material at the same density. The propagation velocity of the shock wave (Uf) through the foam 
was measured and compared with experimental results. The simulations show approximately the 
same agreement with experimental results for U / and agreed with estimates from different equations 
of state and simulations using molecular dynamics. Behavior of the foam ahead of the shock wave 
is also discussed where the porous nature of the foam allows for the formation of a vapor precursor. 

PACS numbers: 52.65.-y, 52.35.Py, 52.35.Qz 
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INTRODUCTION 

There is interest in using low-density metallic foams 
at the National Ignition Facility (NIF), but little detailed 
understanding of the behavior of such foams at conditions 
accessible at NIF. It is unknown whether current hydro- 
dynamic simulation codes and equations of state (EOS) 
can accurately describe the behavior of such foams. 

This work explores the behavior of low-density copper 
foams in regimes of pressure, density and temperature (p, 
p, T) accessible at the NIF using simulations performed 
with HYDRA HYDRA is a single fluid, multi-block, 
multi-material ALE radiation hydrodynamics code. The 
radiation transport was treated in the diffusion approxi- 
mation using a single energy group. Opacities were cal- 
culated from the LEOS tables. The equations of state 
(EOS) for all materials used were from a combined ana- 
lytic / Thomas-Fermi EOS which uses a modified Cowan 
model for the ion EOS, and uses a scaled Thomas-Fermi 
table for the electron EOS Q. 

The results from the HYDRA simulations are com- 
pared with experimental data, the results from molecular 
dynamics simulations of the same system | and calcu- 
lations based on an ideal vapor equation of state and a 
Van der Waals equation of state. 



SIMULATION DETAILS 

The system studied was a slab of copper foam with 
a density of po = 1.3407 g/cm 3 (or 15% of the normal 
density of copper pbuik = 8.938 g/cm 3 ) followed by a 
copper void as shown in Fig. [TJ 

The foam slab had the dimensions 36.8 x 10~ 7 cm 
x 36.8 x 10~ 7 cm x 410 x 10~ 7 cm. The foam was com- 
pressed by a piston using a constant velocity source in 
HYDRA that was incident from x — and moving in the 



+x direction. Piston velocities of U p — 0.2, 2.0. 20.0 and 
50.0 km/s were used in the simulations. The copper void 
was copper with a density of 10 -5 g/cm 3 and was used 
to approximate a vacuum in HYDRA. 
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FIG. 1: A schematic of the set up of the simulation is shown 
above. A piston, moving from left to right, compresses the 
low-density copper foam. Downstream of the copper foam is 
uniformly distributed copper at a density of 10 _5 g/cm 3 to 
approximate a vacuum. 

Two sets of HYDRA simulations were performed. One 
set was three-dimensional and the copper foam was com- 
posed of a tangle of standard density copper filaments 
of cylindrical cross section with a diameter of 4.0 x 10~ 7 
cm; approximately 14 copper atoms reached across a fil- 
aments diameter (the radius of a copper atom is rc u ~ 
1.45 x 10 -8 cm). There were approximately ten zones 
across the diameter of a filament, thus there were approx- 
imately 1.4 atoms in each zone of a filament. Between 
the filaments was the copper void. Figure [5] shows a sec- 
tion of the low-density copper foam. This kind of copper 
foam was first studied using molecular dynamics [jj. 

In the directions transverse to the velocity of the pis- 
ton, there were 92 zones with lengths of 4.0 x 10 -8 cm. In 
the foam region, the length of a zone in the longitudinal 
direction of the piston velocity was also 4.0 x 10 -8 cm. 
In the copper void region the zone lengths were larger 
in the longitudinal direction, but the same in the trans- 
verse directions. The simulations had approximately ten 
million zones. 

The second set of HYDRA simulations were performed 
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in two dimensional RZ geometry. The copper foam was 
simulated as a material that was completely uniform at 
density po- The zoning in the direction along the velocity 
of the piston was the same as in the three-dimensional 
simulations. 




FIG. 2: A section of the copper foam. 

The simulations were performed in parallel with 64 to 
256 processors depending on the number of zones in the 
simulation. Simulations with a piston velocity of U p = 
20.0 km/s and foam slab lengths of k oa , m — 41.0 x 10~ 6 
cm are discussed in detail and results for other values of 
Up and Zfoam reported. 



SIMULATION RESULTS 
Front velocity 

The main goal of this work was to determine the front 
velocity Uf that resulted from the compression of the 
foam by a piston moving at a given velocity U p in order to 
test the understanding of the equation of state (EOS) of 
low-density copper foam and compare with experimental 
results. 

Figure [3] shows, schematically, the compression of a 
uniform material and the resulting shock front. Given 
the velocity of the piston U p , the initial pressure, density 
and temperature (po, Po, Tq) of the material and some 
knowledge of the EOS of the material, one can use the 
so-called jump conditions to determine the velocity of the 
shock front Uf and the pressure, density and temperature 
(p c , Pc, T c ) of the compressed material 
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FIG. 3: A schematic of a uniform material with initial con- 
ditions pa, po, To compressed by a piston (thick vertical line) 
moving with velocity U p . The compression produces a shock 
front moving at velocity Uf and compresses material to the 
state pc, pc, Tc. 



The jump conditions arise from consideration of the 
conservation of mass: 



Pc {Uf - U p ) = p U f 
conservation of momentum: 

PaUfUp =p c -po 
and the conservation of energy: 



p Uf [e c - e 




p c Up. 



(1) 



(2) 



(3) 



Here e is the specific internal energy. 

An experiment that measures Uf for a given U p for 
a material where the initial pressure, density and tem- 
perature values are know then provides some test of the 
knowledge of the EOS of the material. To that end, 
both experiments and simulations commonly measure the 
front velocity as a function of piston velocity. 

For the simulations here, the initial density of the foam 
was po = 1.3407 g/cm and the initial temperature was 
T = 2.56915 x 10~ 5 keV (room temperature). 

The front velocity was determined by dividing the 
length of the foam by the time at which the front breaks 
out of the foam region and into the vacuum region £bo 



Uf = 



'foam 



(4) 



The time at which the front breaks in to the foam region 
is t = 0. The break out time was defined as the time at 
which the density in first zone in the vacuum region was 
equal to the uncompressed foam density. 

Values of U p = 0.2, 2.0, 20.0 and 50.0 km/s were used 
in the simulations. The values of Uf for those simulations 
are shown in Fig. |4] and listed in Table Q] Error bars listed 
for the front velocity arise from the width of the zones 
in the simulations and the time intervals at which data 
from the simulations was recorded. 

Figure 4 shows the results of these simulations com- 
pared with experimental results, the Hugoniot for cop- 
per at 100% pbuik and 15% pbuik from the LEOS ta- 
bles, molecular dynamics simulation results as well as 
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FIG. 4: The front velocity is shown as a function of piston 
velocity: Uf (U p ). Solid black triangles show experimental 
results. Solid red circles show HYDRA filament simulation 
results. Solid blue squares show the HYDRA uniform sim- 
ulation results. Solid magenta stars show the results form 
molecular dynamics simulations. The violet curve shows the 
LEOS table data for 15% normal density copper. The cyan 
curve shows the LEOS table data for normal density copper. 
The green, dashed curve shows the ideal vapor approximation 
for copper at 15% normal density. The green curve shows the 
results using an ideal vapor EOS. The black curve shows the 
results using a Van der Waals EOS. 



estimates of Uf (U p ) using other equations of state and 
the jump conditions. For the piston velocities relevant to 
NIF (Uf > 20 km/s), the simulations (for both uniform 
and filament foam) roughly agree with the LEOS table 
for 15% pbuik and with experimental results [4j. 

At piston velocities of less than about 10 km/s the 
difference between copper at 100% pbuik and 15% pbuik 
from the LEOS tables grows quickly. Both HYDRA sim- 
ulations agree more with the 15% pbuik from the LEOS 
tables than with the 100% pbuik from the LEOS tables. 
This is to be expected for the HYDRA uniform foam sim- 
ulation since that simulation was of uniform copper foam 
of density 15% pbuik whose behavior was governed by the 
LEOS tables. 

The agreement between the HYDRA filament foam 
simulation and the 15% pbuik from the LEOS tables is 
more significant. It illustrates that the behavior of a 
more physically realistic foam is well described by the 
LEOS tables in regimes important for NIF experiments. 

The stars in Fig. show the results from the molecular 
dynamics simulations for values of the piston velocity of 
U p = 1,2.5,5.0,10.0,12.5 and 20.0 km/s 3]. Those re- 
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N/A 



TABLE I: Simulation results for the front velocity as a func- 
tion of piston velocity for a copper foam at 15%po- 



suits show the same general trends as the HYDRA sim- 
ulations and the LEOS table for copper at a density of 
15% pbuik- Results for the molecular dynamics simula- 
tions of the foam (p Cl p Cl T c ) are also shown in figures 
00and0 

Also shown in Fig. and listed in Table U are results 
with an ideal vapor (IV) EOS for copper at (pq, po, To) 
used in the jump condition, equations (fl}, © and ([3]). 
The ideal vapor EOS is 



and 



A ' 



3 P 



(5) 



(6) 



Here A is the atomic mass of copper and is taken to be: 
A = 1.055 x 10~ 22 g. (7) 

The results for Uf using an ideal vapor EOS shown in 
Fig. and listed in Table U were determined by solving 
equations ©, ©, ©, ©, © and © using the initial 
conditions of the foam and the values of U p . Results for 
the final state of the foam (p c , p c , T c ) as determined from 
the ideal vapor EOS are shown in figures [7] and 

Finally, shown in Fig. and listed in Table U are results 
with a Van der Waals (VdW) EOS for copper at (p , p , 
T ) used in the jump condition, equations ((J), © and 
(0). In reduced form, the Van der Waals EOS is 



P = P* 



8 T_ 

3 T, 

E±. _ I 

P 3 



3(£ 



(8) 



and 



e= ttt — -4 1-60210000013 x 10~ 22 . (9) 

2 A 64 p* A 2 w 

Here the critical point of copper is taken from the LEOS 
tables to be: 

T* = 5.27308 x 10" 4 keV, 
p* = 1.98794 g/cm 3 , 
jj* = 5.924738 x 10^ 3 Mbar. (10) 
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The radius of a copper atom according the the Van der 
Waals equation of state and the values of the critical 
point is r-Vdw = 1.62 x 1CT 8 cm. 

The results for Uf using the Van der Waals EOS shown 
in Fig. |4] and listed in Table U were determined by solv- 
ing equations ©, ©, ©, ©, ©, ® and (HU]) using 
the initial conditions of the foam and the values of U p . 
Results for the final state of the foam (p c , p c , T c ) are also 
shown in figures [5j [7] and [8] 

The agreement between the ideal vapor, Van der Waals 
approximations, the HYDRA and molecular dynamics 
simulations and LEOS tables indicates that any reason- 
able estimate of the EOS will lead to a similar U p (Uf) 
dependence. 

Figure |4] and Table UU also show experimental results 
The initial densities of the foam in the experiments 
are not a perfect match for the foams simulated by HY- 
DRA and molecular dynamics. However, the experimen- 
tal results show similar U p (Uf) behavior to that observed 
in the LEOS tables (~ 8% difference) and even the ideal 
vapor approximations (~ 4% difference). The Van der 
Waals approximation shows U p (Uf) behavior that are 
up to ~ 30% different from the experimental results. 
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TABLE II: Experimental and theoretical results for the front 
velocity Uf as a function of piston velocity U p . 



Density 

In the simulations as the piston moved from left to 
right it compressed the foam. In the case of the HY- 
DRA simulations of the uniform foam a shock wave was 
launched in the foam. The foam was compressed to a 
density of p c — 73.9 ± 0.1% pbuik in the compressed re- 
gion at break out time. Here the compressed region is 
defined as the region of the foam between the piston and 
the shock front. Figure [5] shows that that level of com- 
pression occurs shortly after the compression of the foam 
begins and a constant level of compression remains until 
the shock front breaks out into the vacuum region (or 
exits the foam region). 

The results for the filament foam simulations were dif- 
ferent. Nearest to the piston the filament foam is com- 
pressed to a somewhat higher density than the com- 
pressed uniform foam simulations. Away from the piston 
the density of the compressed foam drops with nearly 



FIG. 5: The density profiles at three different times are shown 
for simulations with U v = 20 km/s. The vertical black line 
shows the position of the piston. The blue lines show the 
results of the simulations with the uniform foam. The red 
lines show the results of the simulations with foam made of 
filaments. The green lines show the results for the ideal vapor 
approximation. The black lines shows the results for the Van 
der Waals approximation. The magenta lines show the results 
of the molecular dynamics simulations. The dotted black line 
shows the density of the foam at time t — 0. 



constant slope to the density of the uncompressed foam. 
The plateau in density observed in the uniform foam sim- 
ulation was never observed in the filament foam simula- 
tion. 

The leading edge of the region of compressed foam, 
hereafter referred to as the compaction front, moved 
faster than the shock front of the uniform foam simu- 
lation which caused an earlier break out time for the fila- 
ment foam simulation and thus a faster front velocity U / . 
The average density of the compacted region at break out 
time was (p c ) « 70 ± 10% pbuik- The average density of 
the compressed foam (p c ) was obtained by dividing the 
total mass of copper between the piston and the end of 
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the compressed region at break out time by the volume 
defined by those boundaries. Here the compressed region 
is defined as the region of the foam between the piston 
and the compaction front. The error quoted reflects the 
RMS of the variation in the density in the zones of the 
compressed foam region. 

The density of the compressed foam predicted by the 
ideal vapor EOS is 5.36 g/cm 3 which is lower than the 
HYDRA simulations and, due to the conservation of 
mass, leads to a faster front velocity. 

The density of the compressed foam predicted by the 
Van der Waals EOS is 3.22 g/cm 3 which is lower than 
both the HYDRA simulations the ideal vapor approxi- 
mation and, due to the conservation of mass, leads to a 
faster front velocity. 

The average density of the foam in the compressed re- 
gion in the molecular dynamics simulation shows results 
that are similar to the HYDRA simulations of the fila- 
ment foam although the p (x) slope for the MD simulation 
is less than the HYDRA simulation. 

Pressure 

The HYDRA simulation of the uniform foam reached 
a mass weighted average pressure of p c — 6.77 ± 0.02 
Mbar in the compressed region. Figure [6] shows that for 
the uniform foam simulation the pressure was reached 
shortly after the compression of the foam began and a 
constant level of compression remained until break out 
time. 

The HYDRA simulation of the filament foam reached 
a mass weighted average value of (p c ) — 7.1 ± 0.8 Mbar 
at break out time. Very close to the piston the pressure 
exceeded that attained by the uniform foam simulation. 

The pressure of the compressed foam is 7.15 Mbar for 
the ideal vapor EOS, 9.19 Mbar for the Van der Waals 
EOS and ~7.1 Mbar for the molecular dynamics simu- 
lation. The molecular dynamics simulations (as well as 
the HYDRA simulations of the uniform foam) exhibits a 
nearly constant pressure throughout the compacted re- 
gion. 

p — p Hugoniot 

The measured values of the density and pressure in the 
compressed region can be plotted in pairs and compared 
to the expected Hugoniot. This is shown in Fig. [7] 

In general, the same results observed in Fig. 0] are also 
observed in Fig. [7j However, plotting the data in this 
manner highlights the disagreement that was not easily 
observed in Fig. |4j 

Figure [7| also shows two methods of determining the 
pressure and density for the HYDRA simulations. In 
one method (shown with empty symbols in Fig. [7]) the 
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FIG. 6: The pressure profiles at three different times are 
shown for simulations with U p = 20 km/s. The colors and 
lines are the same as in Fig. [5] 



average values of the pressure and density in the com- 
pressed region were measured directly from the density 
and pressure profiles of Figures [5] and [5] In the other 
method (shown with full symbols in Fig. [7J the density 
and pressure values were determined using the jump con- 
ditions and values of the piston and front velocities (U p 
and Uf). The agreement between these two estimates 
(within error bars) confirms the conservation of mass and 
momentum. 

Generally Fig. [7J shows that the two HYDRA simula- 
tions agree with each other and are close to the LEOS 
tables for copper at a density of 15% pbuik for pressures 
above 0.1 Mbar and are close to the experimental results. 
At the lowest piston velocities the agreement between the 
HYDRA simulation of the uniform foam and the LEOS 
tables is poor. The disagreement (^simulation > Ptabic and 
^simulation > Ptable), at the lowest piston velocities, be- 
tween the LEOS tables and the HYDRA simulations of 
the filament foam may arise due to the vapor precursor 
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FIG. 7: Shown above is the pressure as a function of density 
of the compressed region. The symbols and colors are the 
same as in Fig. [4] Empty symbols show direct measurements 
of the pressure and density while full symbols show estimates 
of the pressure and density using the jump conditions and 
values of the piston and front velocities. 



taking away some of the mass of the copper, thus lower- 
ing the density. The vapor precursor also increased the 
temperature of the filaments (see Fig. [5J and increased 
the pressure. 

The Van der Waals approximation and the molecular 
dynamics simulation show very similar behavior with an 
offset in the density. This offset arises from the differ- 
ences in sizes of the atoms in each system. As stated 
above, the radius of a copper atom in the Van der Waals 
approximation is 1.62 x 10 -8 cm while in the molecular 
dynamics simulation it is 1.45 x 10~ 8 cm. The larger 
atoms in the Van der Waals approximation lead to a 
lower limiting value of the density in the low pressure 
(p 10~ 2 Mbar) region than for the molecular dynamics 
simulation. For the high pressure region (p > 1 Mbar) 
both the Van der Waals and molecular dynamics simula- 
tion tend to constant values of density just as the ideal 
vapor. In that respect, in the high pressure and lower 
pressure region, both the Van der Waals and molecular 
dynamics simulation exhibit ideal behavior. The HY- 
DRA simulations and LEOS tables do not exhibit this 
sort of behavior over the range of pressures and densities 
examined in this work. 
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FIG. 8: The temperature profiles at three different times are 
shown for simulations with U p = 20 km/s. The colors and 
lines are the same as in Fig. [5] Results for the ion and electron 
temperatures from the HYDRA simulations and estimates for 
an ideal plasma are shown. See text for details. 



Temperature and ionization 

In all the HYDRA uniform simulations in the com- 
pressed region the ions were heated to 0.01234 ± 0.00004 
keV and the electrons to of 0.01235 ±0.00005 keV. In the 
HYDRA filament foam, the average values of the tem- 
peratures in the compressed region were 0.0134 ± 0.0007 
keV for the ions and 0.0133 ± 0.0005 keV for the elec- 
trons. All temperatures given above are mass weighted 
average results. 

The values of the temperature for the ideal vapor ap- 
proximation (0.08875 keV), the Van der Waals approxi- 
mation (0.08771 keV) and the molecular dynamics sim- 
ulations (0.077±0.003 keV, the average and RMS of the 
temperature in the compacted region at time t = 15 ns) 
are all similar and nearly an order of magnitude larger 
than the HYDRA results. It is shown below that this 
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is due to the presence of ionization in HYDRA and the 
absence of ionization in the other systems. 

An ideal approximation can be made for a plasma just 
as for a vapor Q. To go from an ideal vapor to an ideal 
plasma, the average number of electrons ionized from an 
atom of the plasma, (Z), is the key. The pressure of the 
ideal plasma depends on the pressure and temperature 
as follows 



Pp 



PP (1 + (Z)) 
A 



(11) 



Here the subscript "P" denotes quantities for the ideal 
plasma. The energy of the ideal plasma is then given as 



E 



-N(l + (Z))Tp + Ne t 



(12) 



N is the number of atoms in the plasma and is the 
energy required to ionize (Z) electrons from an atom. 
Values of e% were determined from the enthalpies associ- 
ated with ionization for copper given in Table Mil [7] . 



Ionization 
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745.5 
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50,200.0 
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61,000.0 


19 


64,702.0 


20 


163,700.0 
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174,100.0 



unionized) vapor and agrees with the behavior shown in 
Fig. [8] 

Using the HYDRA simulations estimates of (Z) (shown 
in Fig. [9]) and the temperature values of the ideal vapor 
(shown in Fig. [8]) gives the results for Tp shown in Fig. [8] 
This estimate of the temperature of the compressed foam 
agrees better with the HYDRA results than do the ideal 
vapor estimates demonstrating that the increased tem- 
perature observed in the ideal vapor approximation, the 
Van der Waals approximation and the molecular dynam- 
ics simulation arises from a lack of ionization. 
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TABLE III: Copper ionization enthalpy values. 

If the ideal vapor and the ideal plasma are at the same 
initial conditions and are compressed by pistons moving 
with the same velocities, conservation of energy allows us 
to combine equations (|6]) and (fl"2j) and solve for Tp, the 
temperature of the ideal plasma as 



Tp = 



T 



(13) 



This shows the the temperature of the (ideal) plasma is 
always less than or equal to the temperature of the (ideal, 



FIG. 9: The ionization profiles at three different times are 
shown for simulations with U p = 20 km/s. The colors and 
lines are the same as in Fig. [5] 



Ahead of the compaction front 

The previous sections described behavior in the com- 
pacted region of the foam. This section describes the 
behavior of the foam ahead of the compaction front. 

The filament foam simulations showed regions of high 
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FIG. 10: Snapshots of a slice through the filament foam simulation at times t = 0, 20, 40 and 60 ps. The top portion of each 
plot shows the ion temperature (keV) while the bottom shows the density (g/cm 3 ); see gray (log) scales at the left. The foam 
is being compressed by a piston moving at 0.2 km/s in the +a;-direction. The filaments are heated well beyond the compacted 
foam region by the vapor precursor. The speed of the vapor precursor is approximately 5.5 km/s, which is a mach number of 
about 1.4. 



pressure ahead of the compaction front. This is evident 
in the middle panel of Fig. |6] which shows the simulation 
at a time of t = 10 ps. At that time the compaction 
front had progressed to a position of nearly x ~ 3 x 10~ 5 
cm, but a region of constant pressure of approximately 
0.5 Mbar is observed between the end of the compaction 
front and about x ~ 3.5 x 10 -5 cm. This region of high 
pressure reflects the extent of the so-called vapor precur- 
sor. The vapor precursor arises because as the piston 
compressed the filament foam, the filaments were heated 
and some copper ablated from the surface of the filaments 
and formed a high temperature vapor. The vapor then 
flowed through the voids in the foam between the fila- 
ments ahead of the compaction front with a velocity of 
approximately 36 km/s for the U p — 20 km/s simulation. 
The vapor precursor heated the filaments ahead of the 
compaction front and increased the average pressure of 
the foam. The uniform foam simulations show no vapor 
precursor since there are no voids for any vapor to travel 
through. 

Beyond the leading edge of the vapor precursor region 
there is another region of constant pressure of about 0.07 
Mbar. The region extends from the end of the vapor 
precursor region to the end of the foam 3.5 x 10~ 5 cm < 
x < 4 x 10~ 5 cm and will be called the preheated region. 
The pressure is elevated in this region because the high 
temperature of the crushed filaments has been conducted 



through the filaments to the end of the foam. All of the 
foam ahead of the compaction front was heated in this 
manner. Examining the temperature profiles shown in 
Fig. |S] shows this behavior as well. The preheat front 
moves through the filament foam at a velocity of approx- 
imately 60 km/s for the simulations with U p = 20 km/s. 
This preheating behavior was observed in only the sim- 
ulation with higher piston velocities. At lower piston 
velocities, the filaments were not heated sufficiently to 
preheat the material ahead of the vapor precursor. 

However, it is not clear that the extent of this preheat 
is entirely physical for the filament foam. The molecu- 
lar dynamics simulations suggest that when heated the 
filaments melt and then the melted portion of the fila- 
ments and are severed from the solid filaments due to 
surface tension effects This behavior arises due to the 
small radius of the filaments (and large surface to vol- 
ume ratio). The melting and severing would preclude 
the heating ahead of the vapor precursor observed in the 
HYDRA simulations of the filament foam. This may be 
expected since HYDRA and the LEOS tables ignore sur- 
face effects. 

For simulations at slower piston velocities the features 
of the vapor precursor and the preheating are combined. 
Figure \TU\ shows snapshots of a slice through a filament 
foam simulation with a piston velocity of 0.2 km/s. Fig- 
ure [11] shows the density, pressure and temperature pro- 
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FIG. 11: From left to right: the density (g/cm 3 ), pressure (Mbar) and temperature (keV) profiles at times t = 0, 20, 40 and 
60 ps. The foam is being compressed by a piston moving at 0.2 km/s in the +a;-direction. The length of the foam slab is 
^foam ~ 23 x 10~ 6 cm. The vapor precursor region is clearly present, however the preheated region is not. The colors and lines 
are the same as in Fig. [5] The vertical dotted line shows the end of the foam slab at time t — 0. 



files for a similar simulation. The snapshots and profiles 
show that the filaments are heated in regions well ahead 
of the compacted foam region (the compacted region for 
these simulation and times is x < 3 x 10~ 6 cm) by the 
vapor precursor. The vapor precursor moves through the 
foam at a velocity of about 5.5 km/s. The filaments then 
melt and flow together forming a liquid with the same av- 
erage density as the foam. This liquid is then compressed 
by the piston. 

Once the vapor precursor exits the foam region it acts 
as a freely expanding vapor. For a U p = 20 km/s the 
front of the vapor precursor in vacuum moves at approx- 
imately 70 km/s. Figure [T2l shows the density, pressure 
and temperature profiles for the vapor precursor in vac- 
uum for Up = 20 km/s. For a U p = 0.2 km/s the front of 
the vapor precursor in vacuum moves at approximately 
2 km/s. Figure [TU1 shows the density, pressure and tem- 
perature profiles for the vapor precursor in vacuum for 



Up = 2 km/s. 

CONCLUSIONS 

The HYDRA simulations, LEOS tables, ideal vapor 
approximation, Van der Waals and molecular dynam- 
ics simulations well reproduce the U p (Uf) dependence 
observed in experiment and HYDRA simulations. This 
indicates a lack of sensitivity in U p (Uf) to simulation 
method or equation of state for copper nano-foams. Sim- 
ilarly, equations of state and methods of simulation have 
little effect on the values of the pressure in the compacted 
region. The density in the compacted region is observed 
to be the same for the HYDRA simulation with filaments 
and the molecular dynamics simulations, but different for 
the HYDRA simulations of the uniform foam and the 
ideal vapor and Van der Waals approximations. Ioniza- 
tion is seen to have a large effect the compressed foam 
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FIG. 12: From top to bottom: the density, pressure and tem- 
perature profiles for the HYDRA filament simulation U p = 20 
km/s at t — 16 ps. The plots show x starting at the t = 
position of the foam/void interface. The rise in all quantities 
after x — 0.7 x 10 -4 cm is due to the vapor hitting the end 
of the simulation at x = 0.75 x 10 -4 cm and bouncing back. 



FIG. 13: From top to bottom: the density, pressure and tem- 
perature profiles for the HYDRA filament simulation U p = 0.2 
km/s at t — 260 ps. The plots show x starting at the t = 
position of the foam/void interface. 



temperature which makes is suitable for testing equa- 
tions of state. The filament foam simulations exhibit a 
vapor precursor that affects the foam and precedes the 
compaction front and is absent in the uniform foam sim- 
ulations. 
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